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1. Objective 

Consider a rigid-link robot with the dynamic model 

f = H(q,p)q 4 C(q,q m ,p)q + G(q;p) 

where *A/*(-) denotes a bounded external disturbance (def- 
initions of other terms and variables can be found in the 
literature)* The objective addressed herein is to find a 
control strategy that exhibits the following features: (1) 
simple to implement , (2) easy to code for program and 
(3) robust to possible time-varying uncertainties 
2- Results 

Define tracking error e a s e = q - q* , where q* denotes 
the desired trajectory. Also let W — t+Dc, y g = q* — Dc, 
and x , = q* - Dc - KK~ l W t where D = D T > 0 and «(t) 
is one of the rate functions which is introduced to adjust 
rate-of-convergen ce (ROC) (see [4]). 

Definition 1 

Let v(t) be defined on [i t 0 , oo). v(t) is in the class Vj if 
i/(t) is nonnegative constant or in V// if v(t) is bounded, 
positive , and decreasing for all t € [^o > 00 )• 

The robust control law is given by 

r= H 9 (q*y)z t +C,{q\q m ;p')y*+ G ,tf 

where K = K T > 0, and #,(.), C,(), and G f (.) are 
simplified versions of H, C, and G , respectively. U a is 
auxiliary control defined by 

Wr? 

Ua ~ \\w\\l + >'{*)’ 

where v(t ) G V/ or V// and r? is a nonnegative scalar 
nonlinear function defined as 

\\H, - H\}\\x,\\ + \\C, -C||||i/,|| + ||<?, -G-M\\ 

< oo||*,|| + ai||«||||»*|| + tt2 = >»- 

In this equation, are constants representing bounds on 
the modelling errors. There are many possible choices for 


t/(<) in U a , and different choices leads to different tracking 
properties. The v(t) defined by 

v(t) = it>i(l + t m Ye~ Vi * n (3) 

are in V/ or V// if m, n, and p are chosen properly (tq 
and t >2 are appropriate positive constants). Specifically, 
if p = -1, m,n = 0, one gets i/(t) = Vie~ v > = pc, 
which gives, U a = ~ \ \ w\\r)+t*i ' is called saturation 
(or boundary layer) controller [1] and has been widely 
used to achieve bounded stability. Another specific choice 
for i/(i) (p = — 1, m = 0, n — 1) is v(t) = V\e 3 , 
which gives the strategy proposed by Dawson, et al. [2], 
77. = . An extreme case, n oo, m — > oc 

and p = -1, gives u(t) = 0, which corresponds to 
the variable structure control [3], U a — “IjWjf 7 ? * s 
shown later, u{t) = 0 gives the fastest convergence, while 
u(t) = po gives the slowest. However, due to physical 
limitations, “too fast” could lead to chattering Hence, 
the choice of u(t) depends on the requirements for ROC. 
transient response, and steady state performance. 

Tracking stability results based on the so-called gener- 
alized energy accumulation principle [4] are given next. 
Theorem 1: Given (1) and (2), if v(t) and «(i) are chosen 
such that 

f i/(r)K 2 (r)dr < C\ < oo Vi € [io.oo), (4) 
Jto 

then stable path tracking is ensured and the rate-of- 
convergence is at least /c -1 (i). 

Proof (outline): The closed-loop model is governed by 
H{q\p)(W + kK- l W) + C(q>q;p)W = - KW 
4 6H(q;p)x s 4 SC{q,q\p)y 9 +6Q(q,t;p) 4 U a . 

where 6H(.) = H t — H,6C( ) = C s — C, 6Q(.) = G s — G — 
X. Introducing the transformation = kW gives 

H(q;p)V 4 C(q,q;p )* = -K * 4 6H(q\p)x 3 K 
4- SC(q,q;p)y,K 4 6G(<ht\p) K 4 U a K. 



According to the criteria in [4], boundedness of the accu- 
mulated generalized energy, f* K^l dr, proves track- 

ing stability. In fact, 


* T KVdr 


, = /■ 
dt 0 

= - / V T H(q-,p)VdT - f V T C{q,q\p)VdT 

Jt a Jt a 

+ f * T {6H(q\p)x,K + 6C(q,q]p)y,K + 6G(q,t-,p)it} 

Jt 0 

+ [ *’ 

Jt o 


\f T U a Kdr. 


The symmetric positive definite property of H ( .) and the 
skew-symmetric property of H ( .) - 2 C(.) yields, 

J c < Cl + f ||*||ipc<fr + / (5) 

where C* = ^HWk 2 |,„,. Inserting (2) into (5), 


J c < Cl+ f ||* \\t)KdT - / * T « 

Jta Jta 

- c*+f\(TW— BS— 

- C V +J ( «T )K imr] + u{T) 

< Cl + f u{r)K 2 {r)dr. 

Jt 0 


Wr, 2 


||W||r? + J/(r) 


dr 


dr 


( 6 ) 


Condition (4) implies J c is bounded. The result follows 

[4]- 

Theorem 2: Given (1) and (2) } if u(t) and «(<) are chosen 
such that 


limsup — - — / z/(r)K 2 (r)dr < < oo Vt € [fo»oo)i 

t — 00 t — to 

then stable path tracking is also ensured . 

Proof: The proof follows the approach used in [4]. 
Corollary If «(<) and u(t ) are chosen such that 

K 2 (t)v{i) < C 2 < 00 , (7) 


then stable path tracking is ensured. 

Proof: Under the condition of the Corollary, it is 

seen that i/(r)« 2 (r)dr < C 2 (< - <o)- Therefore, 

limsup,.^ = iimsup,_ 00 + Cj. 

Several observations are made. First H t1 C $ and G t 
axe not based on q , q and p, but on the desired path 
{q* , q * } and parameters p* which can be precomputed off- 
line. Second, one does not need to re-organize the robotic 
dynamics (so as to isolate unknown parameters) before 
calculating the control torque. Also a simple way to get 
#,(.), C,(.), and G # (.) is to set, H $ = 0, C a = 0, and 
G, = 0, the control torque reduces to r = - KW +U a . 
This gives the same structure as in [2]. However, since 


||£I fl || < 77 , Hi = 0 t C, = 0, and G $ = 0 leads to a larger 
U a which could require more control energy. 

3. Synthesis Examples 

Example 1: (Natural Rate- of- Con verge nee ) 

Assume that a natural ROC is sufficient (*(<)=!). 
Then t 

J c < Cl + f u{r)dr. 

Jta 

Suppose that v(t) is chosen such that J e < C 2 -f- 
J* i/(r)dr < J* , where J* is a design specification. 
If j/(<) = , with v\ > 0,t?2 > 0, then J& ^ 

Cj + 2-Le -Va *°. In order to meet the specification, Vi and 
v? are determined such that C 2 + ^e“ 1 ' 3 * 0 < J* . Suppose 
<o = 0 and the initial condition is such that C 2 = 10. It 
the performance specification is J* = 12, then < 2. So 
by choosing v? > 0 and 0 < v\ < 2^, ^ «!*• 

Example 2: (Variable Structur e Control) 

For any #c(<), is ensured to be less than or equal to 
Cl (see (6)) by choosing v(t) = 0. This implies that the 
ROC can be arbitrarily fast and the accumulated tracking 
error is smaller than any other choice of v. So one might 
conclude that variable structure control gives the best 
control performance and the greatest ROC. However, it is 
its fast speed that causes chattering. So from a practical 
point of view, one should not require too large a ROC 
over the entire period of operation. A piecewise ROC 
may be useful. This can be achieved by methods similar 
to those given in [4], 

4. Comment 

Clearly v plays an interesting role in these robust con- 
trol strategies. First, v is related to the overall tracking 
performance in that the bound on J c depends on the 
choice of u. Second, v specifies the boundary layer in 
the strategies. Since v is time varying, the boundary 
layer is also varying. This property can be used to retain 
the merits of the VSC strategy and avoid the problem of 
chattering. 
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Abstract 

Based on the generalized energy accumulation princi- 
ple, criteria for system stability and performance analysis 
are established in the first part of this work [8]. These 
criteria are of immediate use in many systems. The main 
purpose in this part of the work is to apply these criteria 
to robotic systems. Both adaptive and robust control are 
investigated . 

1. Introduction 

The concept of stability concerning a dynamic sys- 
tem is always important to system engineers. Motivated 
by the generalized energy accumulation principle, cer- 
tain criteria for testing system stability are proposed in 
[8], As a continuation to that work, this paper demon- 
strates the applications erf the established criteria to a 
practical system — robotic system. Specifically, the path 
tracking control problem of robotic systems is consid- 
ered. By introducing the concept of rate transforma- 
tion, new adaptive and robust control strategies are de- 
veloped which achieve stable path tracking and provide 
a priori information about how fast the tracking errors 
will converge. With these strategies, different rates-of- 
convergence (ROCs) can be obtained by simply choos- 
ing a different scalar rate function *(*) (to be defined 
later). It turns out that the conventional adaptive con- 
trol scheme is a special case of the proposed strategies 
(i.e., with K(t) = 1). Global exponential path tracking is 
easily achieved by simply setting n(t) = e Xt , where A > 0. 
Moreover, one may obtain other types of path tracking 
than asymptotic and exponential tracking by choosing a 
proper *. 

2. Review of Stability Criteria 

For convenience and continuity, stability criteria estab- 
lished in [8] are summarized first. 

By introducing a rate function *(<) (see $3 for defini- 
tion), stability criteria represented by integral inequalities 


in [8] can be unified as follows, 

J c = [ K 2 (r)G c [r(r)]dr < C\ < cc, 

Jto 

where G*(.) is a generalized energy function of the sys 
tem and Ca is a real number. To analyze system stabiliu 
and performance, one only needs to verify these equality 
(•ee [8] for more details). It is seen that asymptotic sta 
bility corresponds to «(t) = 1, while exponential stabiliu 
corresponds to k(1) = e Xi . 

3. A Useful Lemma [8] 

Tbe objective is to derive new adaptive and robust con 
trol strategies, which ensure stable path tracking and al- 
low adjustable ROCs. To this end, the concepts of rate 
function and rate transformation are introduced. 
Definition 3.1 Rate Function 

A real function of time, /c(t), is a rate function (denote 
by K(t) 6 S), if it satisfies tbe following conditions: 

(1) n{t) is positive for all t £ [l 0 , oo), 

(2) ic(to) is bounded, 

(3) K(t) is increasing, and 

(4) k(t) is well-defined for t € [to , oo) . 

Note that under these conditions, such a n(t) is invertibh 
and #c _1 (<) is upper bounded and decreasing. Obviou- 
examples for such a rate function include « = 1 , * — I + f 
k = e xt , (1 + t)e u , (l+t)+e A ‘ (A > 0) etc. (see [8] fo 
other types). 

Definition 3.2 Rate Transformation 
Tbe rate transformation is defined as 

* = *(t)x. ( 31 

where *(<) is a rate function as defined before. 

The terminology “rate transformation" is motivated b; 
tbe fact that such a transformation affects the rate o 
convergence of the system, as is shown in the following. 
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Consider a nonlinear system described by 

x = /(x;p). x(<o) = xo, 


Concerning this model, the following structural proper! \ 
iriueAd. ~ i : • 




where x € fZ" i* the system state vector, p € R? » th* a rigid robot with die,, dy^Jmcs (4.1) in which 

system parameter vector, and / € /2" is a nonlinear vac* - C{q, q~, p) is degned as in. [5], then — • ' 

tor function of x and p. Applying the rate transformation : " - v « ** •. ’•* , _ ■ *r- 

(3.1) to (3.2) leads to G(q;p) = ^q,q, V , *)p (4-2) 

♦ = k(t)K- l * + Kf(*f l ;p) i *** - r.,- - -- ' - 

♦(to) = *(<o)x(<o) (3-34) 

where F(k, k, ¥ ,p) is a new nonlinear function. It Aould 
be noted that the transformed system (3.3) is not equiv- 
alent to the original system (3.2) in the sense that the 




stability of (3.2) does not necessarily guarantee the i 
bility of (3.3). However the stability of (3.3) absolutely 
guarantees the stability of (3.2). This is specified by the 
following lemma. 

Lemma 3.3 

Jet k be a rate function as defined before . Consider the 
systems (3.2) and (3.3), related by (3.1). If system (3.3) 
is stable (bounded or asymptotically stable ), the system 
(3.2) converges . to its equilibrium point with a certain 
ROC. . - 

Proof: 

Noting that x = k' 1 * and k" 1 is decreasing, the result 
can be easily obtained. § _ _ . 

Based on this lemma, it is seen that if the traiiffonlfcd 
system (3.3) is at least bounded stable, then the 
one exhibits enhanced stability. This result 
extension of [6] (where a special choice l&r k is utilised, 
i.e., k — e At ) is used in the solution of the tradti&g prob- 
lem. Note that since x could either be ttackipjg VT&ii 
state estimation error, regulation error, or ii 
ing error, the idea behind this lemma could 
in these cases. Investigation of this possibility, beybhd 
the scope of this work, represents an interesting further 
research effort. _/ : ^ 


v 7 [|^U;iO-2C(f f ^p)|v== 0, Vv e R n > (4.3) 

where x € /P*,p € •••>♦«] S 

and € JJ*. It is worth motioning that in (4.2) 
y sb * is nc4 required- This restriction is typically as- 
mined in the literature. Equation (5.3) represents th«- 
well-known skew-symmetric nfeerty. . 2, 

Note, that this equation 

and strqngfy coupled «yst< ^ ^ 

formation concerning the sjfctemperame&rs is availabl 
a priori; Thus the path tfaudcing eontfcJ problem is stale;! 

^riveft the cfatind path, 


___ a fehly nonlinear 
" “ not all the in- 


fcrque, r(f , 

velocities, suchT&it thefic 

f " ' 


Ml 


nd a contro' 
ly on the joini 

path, {q,qj. 

I wjt$ a controllable 
a meters rep - 


— ttyJF J— q 



re jjflC essential difference between this 
topically seen in 
ftatement not only requires 

coavergenceol the tracking errors, but ako imposes the 
requirement d an adj ustab l e RQQ. ybis4s of particular 

41 kdspi with Adjustable ROC 

pork are applies 
ms for robotb 
with control 
defined as 



the tracking error 

4. Application to Robotic Systems _ ,#v; v--. “T-— f’ **• ^ ■ 

The formulation for the dynamics °f serial-link robot . t; fip present an ^pdl iary va^ble defined by 

with n joints is ~ 

H(q,p)q + C(q,q; p)q + G(q-, p) = r, « (4.1) 

where ' 


r € PC 

V>4>4 6 R? 

pel? 

H(q,p)€ir** 

C(q,q,p)q € FT 

G(q;p)€Jr 


- ‘aodltt W present ah amqliary va^le 

\t‘-r : f .\ .r jg) ‘‘•it Vz^T.' " V ‘. uyiPr'v 

- x v&ekr- -m -M* 2w 

f fe- r'S* fit l- SS.’"--.' 

“ i4# r%m^- 

a*d aftVis iii Vactios. Uto 


(4 O 


(4.6; 

(4.7) 


ive definite matrix 
r structure for 


control torque, 

joint positions/velocities/accelerations, 
equivalent system parameters, 
symmetric, poaitive, definite 
matrix, 

joint torques and forces due to 
Coriolis and centrifugal effects, and 

torques ind forces due to gravity. t = H(q\ p)x, + C(q, q;p)y $ + G(q\p) — KW, (4.8) 
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where K — K T > Oiii controller gain matrix, and 
p £ ff ia the vector of parameter estimates as determined 
by the following algorithms. 

Estimation Al gorithm 1 Q-Estimate 1 

pi = -a, J W T *i(q,q,y.,z,)K 2 (T)dT + pi(t 0 ), 

~~ * *'=1,2,...,#, (4.9) 

where a,- > 0 is an adaptation gain (design parameter), 
Pi(t 0 ) i# the initial parameter estimate, and is the i-th 
element of the vector $ defined by Equation (4.2). 

Estimation Algorithm 2 f PI- Estimate } 

Pi = pi(* 0 )-a< / W T $i(q,q t y.,x,)K 1 (T)dT 

Jto 

- PiW T *i(q,q,y t ,z t )K\t), * = 1,2,...,#, (4.10) 

where Pi > 0 is one other adaptation gain. 

Theorem 4.1 

Consider the robot dyn&mics (4.1) with the controller 
(4.8). If the parameters are estimated by (4.9) or (4.10), 
then stable path tracking control is ensured. The rate-of- 
convergence is specified by k~ * (t) , which can be chosen 
by the designer. 

Proof: Combining (4.1) and (4.8) yields the closed-loop 
system dynamics, 

H(q;p){W + kK~ 1 W) + C(q,q;p)W = -KW 

+ 5Z(Pi “ Pi )*<(«> (4 H) 

1=1 

where Property 4.1 was used. Introducing the rate trans- 
formation, * = kW, (4.11) becomes 

9 

H(q\p)*+C(q,q\P ')* = £(P<-Pi )*,(?. 9, *,)t-K 9 . 

1 

According to the criteria in [8], it is only necessary to 
show the boundedness of the accumulated generalized 
energy J e - // * T K^dr . Using the above equation, 
Property 4.1 an3 integration by parts, yields 

J e = - f * T H(q-,p)*dr - / ♦ T C(?,?;p)’i , dr 

Jt, Ju 

,=1 

< [ (pi-Pi)w T *i Kl ( r ) dr 

im 1 ^ 

where Cl = ^HWk 2 | t „,. Inserting (4.9) and apply- 
ing the following relation 

It It, = 5 [/ *^ 7 ^ 7 ] ’ ( 412 ^ 


J c reads 


i= 1 L 

A [Pi ~ PiUo)] 3 

h 2a * 


( 4.13 


1 = 1 


where h = /,‘ # W T 9 i K\r)dr. The boundedness of ./ 
implies that ♦ is at least L 2 . Note that W = th> 

reflult follows. (The result for Estimation Algorithm - 
can be shown in the same way.) § 

4.2 Illustrative Examples 

To make the foregoing concepts clear, three example 
are presented in this section. 

Example 1 Asymptotic Convergence 

Suppose the control torque is of the same structure a> 
in (4.8). If sc = 1, then e, W are defined as before 
and 

*, = - DL 

Estimation Algorithm 1 becomes 


Pi = -ai f W T ii(q,q,y$,x,)dr + p,(<o),» = 1,2, ...,s 

Jto 

Considering the proof of Theorem 4.1, convergence foi 
this case is asymptotic. 

Example 2 

Choosing a(f) as n(t) = 1 + 1 gives 


and 


I, = q- -Di- 


W = ( + D( = 


1 


1 + 1 
* 

1 + t' 


- W . 


Estimation Algorithm 2 is now 


Pi = -cti! W 7 &i(q, q, V», *.)(1 + r) 2 dr 
Jto 

- d.TP T $.(fl,4,V.,l J )(l + 0 2 +P,(<o). 


i = 1,2,...,#. In this case, convergence is stronger thai 
asymptotic due to the choice of k. 

Example 3 Exponential Convergence 

Let k(() be the exponential function *(t) = e A ‘. In thi 

case 

x, = q* - Di - \W 

W = i + D( = '*c- Xt . 



OIRiQiMAL PA OF {$ 

OF POOR QUALITY 


The parameter estimation algorithm is 

pi = -a, / IV 7 $<(?,?, y,,r,)e SAr dr 

Jt, 

- PjW T *i{l, g,y..»«)e 3A< +P<(<o),« = 1,2,...,*- 

This corresponds to the exponential tracking. 

It is observed from the above examples that for differ- 
ent different ROCs for the filtered tracking errors 
are achieved. It is interesting to note that «(*) = 1 cor- 
r«ponds to the conventional adaptive control [l-2][4-5]. 
Note that in this case, W — ¥ and the ROC of W is 
not adjustable. As for the exponential tracking, one only 
needs to choose «(f) = e xt , where X > 0, Also it is possi- 
ble to change the ROC over different time intervals. This 
can be done by the technique shown in [8]. 

4.3 Robust Tracking with Adjustable ROC 

Robust control of robotic systems has been extensively 
investigated recently [1]. Most of the strategies are based 
on upper bounds of the uncertain model. Obtaining such 
bounds, however, is not trivial because H, C, and G are 
complicated matrices depending on q,q and p. Improper 
determination of such bounds may lead to instability. A 
strategy based on the maximum absolute value of each 
element of H, C, and G is suggested as follows. 

Let H t (), C,(.) and (?,(.) represent simplified versions 
of C(.) and G(.), respectively. Also let p* represent 
the nominal system parameters and q 9 and q* represent 
the desired trajectory. 

For the following development, let 

SC = [6dj] = C{q,q\p) - c t {q*,q*;p m ) t 

6G = [Sgi] = G(q\p)-G.{q'\p 9 ). 

The robust control torque is given by 

r= HM*;p*)z t +C,{q\q';pnv,+G,{q m ,p m )-KW+U a} 

(4.14a) 

where K = K T > 0 and U a is an auxiliary control defined 
by 

u a = 

i=t ;=1 

+ (4.144) 

r=iy=i <= i 

In equation (4.14b), U(i,j) are ‘0 — V matrices [9] and 
s#,s§ and sf are scalars to be defined later. The track- 
ing stability of the system is now addressed by the fol- 
lowing result. 

Theorem 4.2 

Let the control strategy be defined by (4.14). If 
B 

9ii \W i *. i \h ij +v(i)‘ 


WiX^c} 


.c _ "•-*> ~Ji 

ij ~ IWy^fo-MO 


«? = - 


M 

\Wi\gi + v(t)' 


where hy = max|fih,j|, ctj = max|5c <; |, y, 
and k end v satisfy 


maxes'; | , 



r?vdr < Ci, < oo 


(4.15) 


or 

k 2 u < Ci, < oo ( 4 - 16 j 

then stable path tracking is achieved. 

Proof (outline): 

Note that with the control (4.14a), the closed- loop 
model becomes 


H{q;p){W + KK-'W) + C(q,q-,p)W = -KW 
+ SH(q\p)x, + 6 C{q,q\p)y, + &G{q\p) + U a . 


Introducing the transformation = kW and using U a in 
(4.14b) gives 




+ C(q,q;p)9 = -K* 

+ EXX + Shij]U(i,j)x,K 

i=ij=i 

+ + 6cij]U(i,j)y,K 

i=i i = i 

+ fy? + 6 9i ]U(i)K. 

1=1 


Considering the performance index J c = f tg *V T f\ 4^/r \- 
is not difficult to show that 


J c 


s: 


1f T K*dT 


< 

+ 

+ 


c; + 


AA [* | W* 

hhX iw**i 


- JvdT 


n n r t 

mi 


l^y <y |c?. 


\w>y.,\cij + v 

gj (* Jr.*™.- K^dr 


+ v 
K 2 udr 


friJt, |Wi| 5 < + *' 


< 

+ 



n 2 vdr 



K 2 vdr + 



K 2 i/dr 


With the choices for k and v as in (4.15) or (4.16), either 
the index is bounded or its time average is bounded and 
the result follows. § 
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This strategy is easy to apply since the upper bound 
for each element, \&hij |m*# > i and can 

be easily obtained by using the facts that jsin(,)j < 1 and 

|cos(.)| < 1. 

Note that there are many possible choices for v (see 
the table below, where /*o > 0, v\ > 0 and vj > 0). The 
impact of u (and k) on system performance is discussed 
later. 

— Table 4.1 Possible Choices for v 


f*o 

Vi 

me-' 1 

»»a+o 

e*2» 

tie - * 3 

«,«-**' 


t-id+i’) 

-uT~ 

1+1 

1+t 



VI .-* »• 




i+< 3 

1+I J 

***** 

Vre-'* 




i+7*~ 

i +< 1_ 

1+t 5 



The strategy presented herein exhibits the following 
features. The structure is simple and most of the re- 
quired computations can be performed prior to real-time 
operation. ,A s for the computation of H $1 C$ and G$ t 
one may choose them to be constant matrices/vector (or 
diagonal matrices for H B and C # ), or simply sero. Addi- 
tionally, time varying uncertainties can be easily handled 
by the strategy. Again since the rate function is utilised, 
the ROC is adjustable. 

4.4 Tracking Performance Analysis 

In addition to the tracking stability, it is important to 
explore the tracking performance that the strategies can 
achieve. The criteria for testing stability can also serve 
this purpose. The following is a brief discussion of this 
issue. Only adaptive control is considered. Referring to 
the proof of Theorem 4.1, it is found that the performance 
index for both the Estimation Algorithms 1 and 2 can be 
computed as 

J C (J, II) = i)o + »?i + fft 

where 

•» = W T H(q(ry,p)WK i {r)\ T=t0 

m = Y ] [ (Pi-Pi)W T *iK\r)dT 

,=l 

m = -H^//(,(r);p)WV(r) | T »«. 

It is seen that only ffr changes for the different estimation 
algorithms. When pi is estimated by the I-Estimate (4.9), 


im t". 
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For the PI-Eatimate, since pi is updated by (4.10), 
Vi = "E / < (P. -w)^ T ^ J (r)dr 

,=l Jt * 

= -E{| J? + lw-M*o)]/<} 

- Pi [ (W T <&) 2 /c 2 dr = r)i(II). 

Jt 0 


Thus 


mill) = min - Pi [\w T ^K 2 (r)dr < 

JtQ 


min- 


Correspondingly it is indicated that 

J C (H) < J C (I) } 

which implies that better tracking performance can be 
achieved by using the Pi-Estimate. This conclusion 
agrees with the comment made in [3]. Simulation results 
presented in Figures 4. 1-4.2 also verify this point (see [7] 
for more details). 

At this point, we are also able to address the effect 
of the initial estimation on tracking performance. Tradi- 
tionally, it is suggested that the initial estimate may be 
-chosen arbitrarily (zero in general). This is because the 
stability is global and the initial estimate does not affect 
tracking stability. However, as clearly shown in (4.13), 
the initial estimate affects the overall tracking perfor- 
mance in the sense that a “better” initial estimate results 
in a tighter bound J c . Simply choosing p, = 0, as sug 
gested typically in the literature, is among the “worst 
choices. Choosing the nominal value of p $ as the initial 
estimate results in a smaller J c > implying better perfor- 
mance. This is also confirmed by simulation results (see 
[7]). These points, however, are not directly evident from 
the Lyapunov stability method. 

Finally, the impact of «(*), and D on system 

performance is discussed. It is noted that K and D are 
required to be symmetric positive definite. Their choices 
are related to the desired robustness, speed of response, 
and disturbance rejection properties. The roles of k and 
v are related to the rate of convergence. Note that since 
the control torque r is defined as in (4.8), if k and v are 
chosen such that e and i rapidly tend to zero, then the 
control torque also rapidly tends to the desired value, r * . 
However, if they are chosen so that the convergence of 
rate is too fast, the control torques in the transition state 
may exceed the admissible values. Hence some trade-offs 
between ROC and control energy have to be made in 
practice. 

5. Concluding Remarks 

This paper has demonstrated the application of the cri- 
teria established in [8] to robotic systems. Performance 
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analysis baaed on these criteria was also given. Addi- 
tional applications of these results can be found in [7]. 
Note that the criteria and their applications are based on 
continuous systems. Given that discrete time systems are 
extensively encountered in practice, extensions of these 
results to discrete-time systems represent an important 
research effort. Due to the limited space, results con- 
cerning this aspect are omitted. Interested readers are 
referred to [7] for details. 

Acknowledgement 

This work was partly supported by NASA grant 
NAGW-2924. 

References 

[1] C. Abdallah, D. Dawson, P. Dorato and M. 
Jamshidi, “Survey of Robust Control for Rigid Robots,” 
IEEE Control Systems Magazine, Vol. 11, No.2, 1991, 
pp. 24-30. 

[2] M. Bodson and S. Sastry, “Adaptive Control: Conver- 
gence and Robustness,” Prentice-Hall, Advanced Refer- 
ence Series, 1989, pp. 229-231. 

[3] S. Gutman, “Uncertain Dynamical Systems — A Lya- 
punov min-max Approach,” IEEE TVans. on Automatic 
Control, Vol .24, 1979, pp. 511-522 

[4] I.D. Landau and R. Horowitz, “Synthesis of Adaptive 
Controllers for Robot Manipulators Using Passive Feed- 
back System Approach,” Int. J. Adaptive Control and 
Signal Processing. Vol. 3, No.l, 1989. 

[5] J.J. Slotine and W. Li, Applied Nonlinear Control, 
Prentice Hall, Englewood Cliffs, 1991. 

[6] Y.D. Song, R.H. Middleton and J.N. Anderson, “Study 
on the Exponential Path Tracking Control of Robot Ma- 
nipulators via Direct Adaptive Methods,” to appear in 
Int. J. of Robotics and Autonomous Systems, Vol. 9, 
1992, pp.271-282. 

[7] Y. D. Song and J. N. Anderson, “The Principle of 
Generalized Energy Accumulation and Its Application to 
the Motion Control of Robotic Systems,” Tech Report 
No. MCTR-0792-02, July 1992. 

[8] Y. D. Song and J. N. Anderson, “System Stability 
and Performance Analysis Based on Generalized Energy 
Accumulation: Part I — Criteria Development,” IEEE 
1st Int. Conf. on Control Applications, September, 1992. 


(a) j«n< I (fc) kssi J 



u-* (•) (•) 


{c)jQiMl 


{4)fafci2 



4.1 Tracking Errors with Adaptive Control 



4.2 Performance Index with Adaptive Control 





